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There is a worldwide effort toward the development of a large TPC (Time Projection Chamber) 
devoted to directional Dark Matter detection. All current projects are being designed to fulfill a 
unique goal : identifying weakly interacting massive particle (WIMP) as such by taking advantage 
of the expected direction dependence of WIMP-induced events toward the constellation Cygnus. 
However such proof of discovery requires a careful statistical data treatment. In this paper, the 
discovery potential of forthcoming directional detectors is adressed by using a frequentist approach 
based on the profile likelihood ratio test statistic. This allows us to estimate the expected significance 
of a Dark Matter detection taking into account astrophysical and experimental uncertainties. We 
show that the energy threshold and the background contamination are key experimental issues for 
directional detection, while angular resolution and sense recognition efficiency only mildly affect the 
sensitivity and the energy resolution is unimportant. This way, we found that a 30 kg.year CF4 
directional experiment could reach a 3<r sensitivity at 90% C.L. down to 10~ 5 pb and 3.10 -4 pb for 
the WIMP-proton axial cross section in the most optimistic and pessimistic detector performance 
case respectively. 



PACS numbers: 95.35.+d, 14.80.-j 
I. INTRODUCTION 



An ever increasing body of evidence supports the ex- 
istence of Cold Dark Matter (CDM) as the major con- 
tribution to the matter budget of the Universe. On the 
largest scale, cosmological measurements [l| tightly con- 
strain the CDM relic density whereas on a local scale, the 
discrepancy between the rotation curves of spiral galax- 
ies, as measured by 21 cm observations, and the ones that 
would be inferred from the luminous matter, indicates 
that spiral galaxies, including the Milky Way, should be 
embedded in a Dark Matter halo 0, Q . 
Directional detection of galactic Dark Matter has been 
proposed as a powerful tool to identify genuine WIMP 
events as such jj). Indeed, one expects a clear asymme- 
try in the angular distribution of WIMP-induced events 
in the direction of motion of the Solar system, which 
happens to be roughly in the direction of the constella- 
tion Cygnus. As the background distribution is expected 
to be isotropic in the galactic restframe, one expects a 
clear and unambiguous difference between the WIMP- 
induced signal and the background one. Recent studies 
have shown that, within the framework of dedicated sta- 
tistical data analysis, a low exposure directional detector 
could lead either to a high significance discovery of galac- 
tic Dark Matter 0-0] or to a competitive exclusion Q. 
Non-directional direct detection of Dark Matter is en- 
tering in a new era as several detectors are starting to 
exclude the upper part of the predicted supersymmetric 
parameter space, either in the spin-independent channel 
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9 12] or in the spin-dependent one (l3rl20j . On the other 
hand, forthcoming directional detection projects are de- 
signed to claim a proof of discovery. Both approaches 
require a dedicated statistical treatment, taking into ac- 
count systematics related to detector performance as well 
as astrophysical uncertainties. Following [IHIIIJ], we use 
a profile likelihood ratio test statistic to assess the dis- 
covery potential of directional detection. 
We give a complete overview of the effect of the main 
experimental issues on the discovery potential of direc- 
tional detection of Dark Matter. Our aim is to establish a 
weighted wish list that could be used when trying to opti- 
mize the design of a Dark Matter directional detection ex- 
periment. The experimental issues considered hereafter 
are: the background contamination, the energy thresh- 
old, the sense recognition efficiency, the angular resolu- 
tion, the energy resolution and the background modeling. 
With a profile likelihood method, we estimate the sensi- 
tivity of a given directional experiment by taking into ac- 
count the most relevant astrophysical uncertainties, e.g. 
the local Dark Matter density, the WIMP velocity distri- 
bution and the velocity of the Solar System orbit. 
The paper is organised as follows. Section |n] presents a 
brief introduction of directional Dark Matter detection, 
both from an experimental point of view and from a the- 
oretical one. Section Mil introduces the definition of the 
likelihood function we have used, as well as the formal- 
ism of a frequentist approach to estimate the significance 
with a profile likelihood ratio test statistic. Section IIVI 
focuses on the application of the profile likelihood test to 
estimate the impact of the different experimental issues 
on the sensitivity of upcoming directional detector. Fi- 
nally, in section [V] we conclude this study by considering 
two very different detectors with their respective perfor- 
mance to evaluate both their sensitivity to Dark Matter 
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from supersymmetric models and their competitivity in 
comparison to existing limits. 



II. DIRECTIONAL DETECTION FRAMEWORK 

A. Detector configuration 

There is a worldwide effort toward the development 
of a large TPC (Time Projection Chamber) devoted to 
directional detection Q and all current projects [25l — 
f30| | face common experimental challenges and share a 
unique goal : the simultaneous measurement of the en- 
ergy (E r ) and the 3D track (tt r ) of low energy recoils, 
thus allowing to evaluate the double-differential spectrum 
d 2 R/dE r dn r down to the energy threshold. This can be 
achieved with low pressure gaseous detectors (TPC) and 
several gases have been suggested : CF 4 , 3 He, C4H10 
or CS2 • As a matter of fact, exp erimental key issues 
include : sense recognition [3ll433| . angular and energy 
resolutions, energy threshold as well as the residual back- 
ground contamination. 

Bearing in mind the need for a directional detector opti- 
mization d, [34| at the current stage of detector design, 
we present in section IIVI a study of the effect of detector 
configurations in the case of a working example : a low 
exposure (30 kg.year) CF 4 TPC, operated at low pres- 
sure and allowing 3D track reconstruction as proposed 
by the MIMAC collaboration [30|. 



The directional recoil rate is given by |3£ 
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with m x the WIMP mass, m r the WIMP-nucleus reduced 
mass, po the local Dark Matter density, <tq the WIMP- 
nucleus elastic scattering cross section, F(E r ) the form 
factor (using the axial expression from 40]), v m i D the 
minimal WIMP velocity required to produce a nuclear re- 
coil of energy E r and q the direction of the recoil momen- 
tum. Finally, /(v m in,<?) is the three-dimensional Radon 
transform of the WIMP velocity distribution f(v). Us- 
ing the Fourier slice theorem [33], the Radon transform 
of the multivariate Gaussian is, 
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As outlined in 043 > a clear and unambiguous signa- 
ture in favor of a Dark Matter detection is given by the 
fact that the recovered main recoil direction is shown to 
be pointing toward Cygnus within a few degrees. How- 
ever, to assess the significance of the discovery, in a fre- 
quentist approach and taking into account astrophysical 
uncertainties, a complete statistical analysis is required, 
as outlined below. 



C. Astrophysical uncertainties 



B. Directional detection 

Directional detection depends crucially on the local 
WIMP velocity distribution [35l - [37j . The isothermal 
sphere halo model is often considered but it is worth go- 
ing beyond this standard paradigm when trying to ac- 
count for all astrophysical uncertainties. The multivari- 
ate Gaussian WIMP velocity distribution corresponds to 
the generalization of the standard isothermal sphere [38[ 
with a density profile p(r) cx 1/r 2 , leading to a smooth 
WIMP velocity distribution, a flat rotation curve and 
no substructure. The WIMP velocity distribution in the 
laboratory frame is given by, 
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where er v = di&g[a x , a y , a z ] is the velocity dispersion ten- 
sor assumed to be diagonal in the Galactic rest frame (x, 
y, z) and vq is the Sun's velocity vector with respect to 
the Galactic rest frame. When neglecting the Sun pecu- 
liar velocity and the Earth orbital velocity about the Sun, 
Vq corresponds to the detector velocity in the Galactic 
rest frame and is taken to be Vq = 220 km.s" 1 along 
the y axis pointing toward the Cygnus constellation at 
(i Q = 90°, b Q = 0°). 



The effect of astrophysical parameters on exclusion 
limits deduced from direct detection, both from energy, 
annual modulation and direction measurements, has 
been investigated in detail in jl, |H, H3, lill-li^]. For a 
complete discussion on their evaluations, uncertainties 
and effect on dark matter detection, we refer the reader 
to [U Ell and references therein. For dark matter 
search, both direct and directional, it is indeed of great 
interest to account for uncertainties on astrophysical 
parameters, as in (23|, although it has been done for the 
escape velocity only, within the framework of a profile 
likelihood ratio test statistic. 

Three astrophysical parameters play a key role : the 
local dark matter density po, the Sun's velocity vector 
Vq and the WIMP velocity distribution, parametrized 
by the velocity dispersions along the three axes a x . y . z . 

As the WIMP event rate is proportional to the quan- 
tity po x Up, the value of the local dark matter density 
po directly affects the estimation of the WIMP-proton 
cross section a p . A standard value po = 0.3 GeV.c~ 2 
is usually used for the sake of comparison of various 
direct detector results [Ii|[50|. However, there are still 
debates on the value of po and several papers, using 
different techniques, have found the following values of 
the local Dark Matter density : po(GeV /c 2 / cm 3 ) equal 
to 0.32 ± 0.07 SI, 0.43 ± 0.11 ± 0.10 M, 0.3 ± 0.1 
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and 0.39 ± 0.03 [53|. Clearly, systematic errors 
arise from uncertainties in modeling the Milky Way. 
Interestingly, S. Garbari et al. '54} point out that it is 
vital to measure the vertical dispersion profile of the 
tracers to recover an unbiased estimate of po. Depending 
on the choice of the velocity distribution of star tracer 
population (from isothermal to non-isothermal), their 
estimation of po can vary by a factor of 10 or so. 
Recently, M. C. Smith et al. [55| have computed a local 
dark matter density of 0.57 GeV/c 2 /cm 3 using SDSS 
data [5(| by studying the kinematics of Galactic disk 
stars in the solar neighborhood. 

Valuable information are also extracted from recent 
high resolution N-body simulations of M ilky Way like 
objects. In particular, it is shown in [57] that the 
estimation of the local Dark Matter density, on the 
stellar disk and at ~8 kpc from the Galactic center, is 
found to be always larger, by 20% or so, than the average 
density in a spherical shell of same radius, i.e the above 
quoted quantity inferred from dynamical measurements. 
Interestingly, the local Dark Matter density distribution 
is predicted [EH to be remarkably smooth, varying by 
less than 15% at the 99.9% C.L. from the average value 
over an ellipsoidal shell. This highlights the fact that 
the local density, at the Sun's location, should not differ 
from the above estimations. 

These arguments, both from observations and simula- 
tions, add still more weights to the need of taking a 
possible range of values for the local Dark Matter density 
rather than a fixed standard value (po = 0.3 GeV.c -2 ), 
although it allows a fair comparison of various dark 
matter search results. Hence, to account for all estima- 
tions and uncertainties, the local Dark Matter density 
Po is treated hereafter as a nuisance parameter in the 
following, considering the range 0.3 ± 0.1 GeV/c 2 /cm 3 , 
see tab. U 

On the choice of mean value and uncertainty, it is worth 
emphasizing that the local Dark Matter density is indeed 
a key experimental issue for directional Dark Matter 
search (as for direction-insensitive one) but its effect is 
rather simple to handle. For instance, taking a larger 
value (e.g. as proposed by f55j). say 0.4 instead of 0.3, 
will mainly change the sensitivity by 33%, giving in fact 
a sensitivity to smaller cross sections. For a study aiming 
at showing the reach of directional detection, a modifica- 
tion of the range chosen for the local Dark Matter density 
does not affect the result. This is not the case however, 
when presenting an experimental result as oulined above. 

The second astrophysical parameter to be carefully 
handled is the Sun's velocity v@ equal to the local cir- 
cular speed when neglecting the Sun's pec uliar velocity. 
The standard value is 220 ± 20 km.s" 1 [59[. As outlined 
in [H, H3, Hil, recent determinations of its value span 
on a wide range and the impact on directionality has 
been studied [1, |fi3]. Several papers, using different 
techniques, have studied the value of local circular speed, 
giving slightly different values. They found Vq (km.s -1 ) 



Nuisance parameters 


Gaussian parametrization 


„„ [HpV lr 2 /rm 3 l 






220 + SO 


a x [km/s] 


220/^/2 ± 20 


G y [km/s] 


220/^/2 ± 20 


cr z [km/s] 


220/^/2 ± 20 



TABLE I: Gaussian parametrization (mean and standard de- 
viation) of the different astrophysical nuisance parameters. 



equal to 254 ± 16 (6TJ, 200 - 280 [62j], 236 ± 11 [H, 
221 ± 18 [El and 218 ± 7 [EE|j. It is worth noticing that 
the evaluation of the Sun's circular velocity depends 
crucially on the value of the distance of the Sun to 
the galactic centre (Rq) which is also poorly known : 
R = 8.4 ±0.4 kpc @ andi? = 8.33 ± 0.35 kpc [g2|. In 
the following, the local circular speed is treated as a nui- 
sance parameter, considering the range 220 ± 30 km.s -1 , 
see tab. Q] 

Finally, uncertainties in the local WIMP velocity dis- 
tribution must be accounted for. The effect of halo mod- 
eling on exclusion limits and allowed regions has been 
studied 0,Hg|- Indeed, recent results from N-body sim- 
ulations are in favor of triaxial Dark Matter haloes with 
anisotropic velocity distributions and potentially con- 
taining substructure as sub- haloes (clumps) and dark 
disk [4ll I69h72i ]. Moreover, recent observations of Sagit- 
tarius stellar tidal stream have shown evidence for a tri- 
axial Milky Way Dark Matter halo [73| . However, it is 
noteworthy that this result holds true at large radius (60 
kpc) and N-body simulations have shown that there can 
be significant variations of the axis ratios with radius 
[74j . S. H. Hansen and B. Moore identified a universal 
relation between the radial density slope and the veloc- 
ity anistropy leading to ft ~ 0.1 at the solar neighbour- 
hood [75|. Using a sample of 1700 solar neighbourhood 
halo subdwarfs from the Sloan Digital Sky Survey [56| . 
M. C. Smith et al. have evaluated the halo velocity dis- 
persion and found the anisotropy parameter j3 to be ~ 0.5 
[76j |. Interestingly, the stellar halo exhibits no net rota- 
tion. 

Even if it should be noticed that halo stars have a pri- 
ori a different density profile from the dark matter, these 
facts can be taken as hints in favor of an anisotropic Dark 
Matter velocity distribution. To account for this effect, 
we consider a multivariate Gaussian WIMP velocity dis- 
tribution, i.e. anisotropic, but without substructures. 
It is parametrized by the velocity dispersions o~ XiV , z . Ef- 
fect of non-smooth halo model with substructures and/or 
streams will be adressed in a forthcoming paper. 
The velocity anisotropy fi(r), is defined as [77j . 
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According to N-body simulations with or without 

baryons [4ll l43l I7a - l80j . the (3 parameter at a radius 
R & ~ 8 kpc from the Galactic center spans the range 
— 0.4 which is in favor of radial anisotropy. Indeed, such 
radial anisotropy is expected as the gravitational poten- 
tial is mainly a function of the radius. In the following, 
the velocity dispersions <J x ,y,z are treated as nuisance 
parameters, considering the range 220/y/Z ± 20 km/s, 
see tab. U which corresponds to (3 = ± 0.25. We 
used a null mean value of the velocity anisotropy for 
the sake of comparison with other experiments, but the 
large uncertainty considered enables us to account for 
anisotropic halo model when computing the discovery 
potential of directional experiments. Moreover, as 
shown in [7|, the effect of an extremely anisotropic 
halo model with j3 = 0.4 can be handled with such 
parametrization of the WIMP velocity distribution, 
as long as it can be approximated by a multivariate 
Gaussian, without introducing bias in the estimation of 
the WIMP properties (mass and cross section). 

In this study, we have not used the escape velocity v csc 
as an astrophysical nuisance parameter, since directional 
detection with low energy threshold is almost insensitive 
to this parameter. Indeed, in the case of a WIMP mass 
of 100 GeV/c 2 the minimal speed to produce a 5 keV 
Fluorine recoil is 130 km.s -1 , far below the median 
value of [8l| equal to 544 km/s. Thus, the impact of v csc 
is relevant for heavy targets and low WIMP mass in the 
so-called threshold region, where the experiment is only 
sensitive to the tail of the WIMP velocity distribution. 
In the following, we have considered an escape velocity 
taken as infinity to simplify the calculations. 

As a conclusion, the use of a profile likelihood ratio 
test statistic allows us to account for uncertainties on 
the astrophysical parameters, as they are treated as nui- 
sance parameters. This approach is a step beyond the 
"standard Dark Matter halo paradigm", i.e. isotropic 
isothermal Dark Matter halo with fixed value of density. 
Evaluating the properties of the Dark Matter halo is in- 
deed still a subject of debates and directional detection 
could bring valuable information, such as an evaluation 
of the anisotropy parameter, as shown in Q- 



in a benchmark case. 



A. The likelihood function 

When considering directional detection of Dark Mat- 
ter, the event distribution is given by the double dif- 
ferential spectrum (eq. [5]) which is a function of both 
the recoil energy and direction. On the other hand, the 
background event distribution is expected to be isotropic 
and will be considered as flat in energy unless otherwise 
stated. In this study, we are also taking into account five 
different astrophysical uncertainties as: the local Dark 
Matter density po, the circular velocity of the Sun vq 
and the velocity dispersion along the three axes of the 
halo {a x , <jy, a z }. Then, considering the extended def- 
inition of the likelihood function (see Ref.[22|]) to take 
into account the Poisson statistics of the number of ob- 
served events N and the likelihood terms associated with 
each nuisance parameters, the full likelihood function is 
defined as, 



N r 
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Ms + Mb 



(5) 



where v = {po,VQ,a x ,a y ,a z } represents the set of nui- 
sance parameters, /if, = Rb x £ and Ms correspond to the 
number of expected background and WIMP events re- 
spectively, where £ corresponds to the exposure. N is 
the total number of observed events, R n refers to the en- 
ergy and direction of each event while the functions S (see 
eqJ3]) and B (isotropic distribution for the angular part 
and following ea !15l for the energy part) are the direc- 
tional event rate d 2 R/dEndQ,R of the WIMP and back- 
ground events respectively. The functions J£i(yi) refer to 
the likelihood function associated with each astrophysi- 
cal nuisance parameter. They are considered as Gaussian 
functions with the parametrization from Table HI 



III. THE PROFILE LIKELIHOOD RATIO TEST 
STATISTIC 

We are interested in estimating the expected signifi- 
cance of a discovery of Dark Matter with directional de- 
tection. In this section, we first introduce the definition 
of the likelihood function that we have used. In particu- 
lar, astrophysical uncertainties are accounted for. Then, 
the formalism of a frequcntist approach to estimate the 
significance using a profile likelihood ratio test statistic 
is presented. Finally, the impact of astrophysical uncer- 
tainties in the estimation of the significance is presented 



B. The profile likelihood ratio 

The significance of a new process, in a frequentist ap- 
proach, is commonly estimated by using the profile like- 
lihood ratio test. We present hereafter a brief overview 
of the subject by recalling basic relations for the reader's 
convenience. We refer to [2l[ for a comprehensive dis- 
cussion. The profile likelihood ratio test corresponds to 
a hypothesis test against the null hypothesis Ho (back- 
ground only) against the alternative Hi which includes 
both background and signal. The strength of the pro- 
file likelihood ratio relies in the fact that one can add 
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some nuisance parameters both from an experimental 
side (efficiencies, resolutions,...) and from a theoretical 
side (signal or background modeling). In our case, the 
Hq hypothesis corresponds to the a p — case (back- 
ground only), while the alternative Hi corresponds to 
the case where a p is taken different from (background 
and signal) . Then, in the case of a discovery, we test the 
background only hypothesis on the observed data and try 
to reject it using the following ratio 



A(0) 



J?(ap = 0, Rb, v) 
J*?(d p ,R b ,v) 



(6) 



As discussed in [21( the test statistic in such case is de- 
fined as, 



qo = 



-21nA(0) d p > 







dp < 



(7) 



As it can be seen from the previous equations, < 
A(0) < 1 and q > 0. Large value of qo implies a large 
discrepancy between the two hypotheses, which is in fa- 
vor of a discovery. The p- value po is defined as 



Po 



f(qo\H )dq Q 



(8) 



where f{qo\Ho) is the probability density function of qo 
under the background only hypothesis Hq. Then, pg cor- 
responds to the probability to have a discrepancy, be- 
tween H and Hi, larger or equal to the observed one 
q ohs . As an example, a p-value of p n = 0.00135 corre- 
sponds to a 3cr signal observation. The main issue here is 
to have a correct estimation of the f{qo\Ho) distribution 
in order to compute the significance for a given data set. 
However, following Wilk's theorem, qo asymptotically fol- 
lows a x 2 distribution with one degree of freedom, see. 
PH - Then, in such case, the discovery significance Z is 
simply defined as Z = y/ q ohs , in units of a, i.e. Z = 1 
corresponds to a significance of 68%. Figure Q] presents 
the probability density functions f(qo\Ho) of qo under 
the background only hypothesis Ho with/ without astro- 
physical uncertainties. It is worth noticing that f(qo\Ho) 
is very well fitted by the \\ distribution. Indeed, in both 
cases: with and without astrophysical uncertainties, cor- 
responding to the blue and black histograms, we found 
X 2 /d.o.f equal to 95.27/99 and 99.79/99 respectively. 



C. Impact of astrophysical uncertainties 

We aim at evaluating the discovery potential of upcom- 
ing directional detection experiment as a function of their 
experimental performance (resolutions, threshold, sense 
recognition) taking into account astrophysical uncertain- 
ties (see sec. Ill Cl) . In the following, we estimate the im- 
pact of the latter by studying two different cases: with 




FIG. 1: Probability density functions of qo under the 
background only hypothesis Ho with/without (blue solid 
line/black dotted line) considering astrophysical uncertain- 
ties, estimated with 10,000 Monte Carlo simulations. The 
red long dashed line correspond to the xi distribution. 



and without uncertainties from astrophysics. Not consid- 
ering the latter corresponds to remove the functions 
from eq. ([5]) and to consider all the astrophysical param- 
eters as perfectly known (without error bars). To do the 
comparison, we ran 10,000 Monte Calo simulations of a 
benchmark model with 100 WIMP events and 50 back- 
ground events with a WIMP mass of 50 GeV/c 2 . The 
distributions of the significance Z corresponding to the 
two scenarios are presented on figure^ The black dashed 
histogram corresponds to the case without uncertainties 
and the blue solid one corresponds to the case where 
the uncertainties are taken into account. One can easily 
deduce that the mean value of the significance E(Z) is 
almost the same, i.e. 9.3 and 9.2 respectively. However, 
as shown on figure [2l the spread of the Z distribution is 
much wider when taking astrophysical uncertainties into 
account. This leads to a large difference between the two 
cases when considering a confidence level on the value of 
Z. Indeed, we defined Zgo as the value of the significance 
that could be at least reached by an experiment 90% of 
the time. Then Zgo is found by solving the following 
equation, 



Zqq 



f(Z)dZ = 0.9 



(9) 



where f(Z) corresponds to the normalized distribution of 
the significance Z. Then, the value of Zgo with/without 
taking into account the astrophysical uncertainties is 5.9 
and 7.9 respectively, leading to a large effect on the ex- 
pected sensitivity to Dark Matter of a directional detec- 
tor. 

As a conclusion, this highlights the fact that it is nec- 
essary to take into account astrophysical uncertainties 
when estimating the sensitivity of a given directional de- 
tection experiment to Dark Matter. In the following, 
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Significance Z 



FIG. 2: Normalized distributions of the significance f(Z) 
with/without (blue solid line/black dotted line) considering 
astrophysical uncertainties. Distributions were generated us- 
ing 10,000 Monte Carlo simulations with 100 WIMP events 
and 50 background events with a WIMP mass m x = 50 
GeV/c 2 . We found mean values of the significance of 9.2 
and 9.3 and a 90% confidence lower limit Zgo of 7.9 and 5.9 
respectively. 
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FIG. 3: Evolution of the mean significance E(Z) as a func- 
tion of the expected WIMP fraction X = N S /(N S + N b ) for 
three different values of the expected total number of events 
Ntot = 100 (black solid line), 50 (red dotted line) and 25 (blue 
long dashed line) . This study has been done by considering a 
WIMP mass of 50 GeV/c 2 . 



we define the sensitivity of a given directional experi- 
ment by the lower bound of the 3c discovery region at 
90% Confidence Level, given by the Zgo = 3 limit in the 
(m x , log 10 (<7 P )) plane. 



IV. OPTIMIZING DIRECTIONAL DETECTION 

All current directional projects [25T - [3Ci| face common 
challenges. As a matter of fact, experimental key issues 
include : sense recognition pul - l33j . angular and energy 
resolutions, energy threshold as well as the residual back- 
ground contamination. 

In this section, we evaluate the evolution of the signifi- 
cance of a discovery as a function of the detector perfor- 
mance. In a first place, we will be interested in the effect 
of the number of WIMP and background events in the 
data. Then, in a second part, we will consider some ex- 
perimental issues, like the energy threshold E t h, the an- 
gular resolution <r 7 , the sense recognition efficiency, the 
energy resolution <7e,. and the energy background mod- 
eling. In the following, we consider the mean significance 
E(Z) and the sensitivity, estimated from the limit corre- 
sponding to Zgo = 3 in the (m x , \og w (a p )) plane, using 
1,000 Monte Carlo simulations for each given set of in- 
put parameters. Unless otherwise stated, we consider a 
10 kg CF4 directional detector with a recoil energy range 
of [5, 50] keV with perfect angular and energy resolu- 
tions and with sense recognition. To study one by one 
the impact of a given experimental issue on the sensi- 
tivity, this parameter will be degraded keeping all other 
ones unchanged. 



A. Residual background contamination 

Zero background is often referred to as the ultimate 
goal for the next generation of direct detection experi- 
ments in deep underground laboratories. However, owing 
to the large intrinsic difference between the WIMP- 
induced and background-induced spectra, directional 
detection could accommodate to a sizeable background 
contamination. In this section, we investigate the effect 
of a residual background contamination, i.e. in the final 
dataset after all selections based e.g. on track-length 
versus energy and fiducialization. Hence, we try to 
account for the effect of an irreducible background 
contamination, which can only be treated by means of 
statistical analysis. 

To begin with, on figure [3] we studied the effect of 
background/signal contribution to data. Figure [3] repre- 
sents the evolution of the mean significance E(Z) as a 
function of A = N s / (N s + Nb), which corresponds to the 
fraction of WIMP events contained in the data, for three 
different total numbers of events N to t — 25, 50, 100. 
Hence, A — > means that the data are background 
dominated. Then, from figure one can observe that 
for any total number of events, the mean significance 
E(Z) increases almost linearly with the data purity (A). 
Also, for a given value of A, increasing the total number 
of events will enhance the significance. In other words, 
for a given value of A = 0.5 for example, the significance 
of a Dark Matter detection can be improved from 3cr to 
6(7 by having an exposure 4 times larger. 

Figure 0] presents the lower bound of the 3<r discovery 
region at 90% C.L. (Zg — 3 limit) in the (m x , log 10 (<7 p )) 
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plane. Two cases are presented : background-free and 
high background contamination (10 kg -1 . year -1 ). For 
convenience, the curves of iso-number of WIMP events 
are presented (dashed lines) in the case of 5, 30 and 150 
WIMP events. One can first notice that a 30 kg. year 
CF4 directional detector may allow to achieve, in the 
background-free case, a 3a discovery of Dark Matter 
down to - 10~ 5 pb at low WIMP mass (~ 10 GeV/c 2 ) 
and down to ~ 10~ 4 pb for high WIMP mass (~ 
1 TeV/c 2 ). Adding a large fraction of background in the 
data (10 kg _1 .year~ 1 ) results in a loss of about one order 
of magnitude in the lower bound of the 3cr discovery re- 
gion. This highlights the fact that directional detection 
can accommodate to a sizeable background contamina- 
tion, noticing that such a background event rate in the 
final dataset is very large. Interestingly, one may com- 
pare the 3cr lower bound to the curves of iso-number of 
WIMP events. At low WIMP mass (~ 10 GeV/c 2 ) about 
5 WIMP events are enough to claim a 3cr discovery, in 
the background-free case, while about 30 WIMP events 
are needed for highly-background contaminated data. At 
high WIMP mass (~ 1 TeV/c 2 ) these numbers reach 30 
and 150 respectively. 

To conclude on this first study, aiming at evaluating the 
effect of N s and Nj, on the mean significance, one can 
easily appreciate the fact that using directional detec- 
tion, strong evidence in favor of a Dark Matter detec- 
tion could be reached even with low statistics and rather 
large background contamination. This low sensitivity to 
background contamination and the possibility to evalu- 
ate the Dark Matter and background components in the 
final dataset, give a major interest for directional detec- 
tion, especially at the present stage when non-directional 
experiments start to observe candidate events whose ori- 
gin is difficult to assess [9l4l^|. On the other hand, we 
emphasize that, for directional detection, an unambigu- 
ous proof of discovery of Galactic Dark Matter would be 
given by a main recoil direction pointing toward the con- 
stellation Cygnus within few degrees, combined with a 
high significance. However, we caution that the presence 
of a stream or another Dark Matter substructure in the 
vicinity of the Solar System could deviate the recover- 
ing angle from the direction of the constellation Cygnus 
of a few de gree s, depending on the relative densities, as 
outlined in 18211. 



B. Effect of the energy threshold 

As for direction-insensitive direct detection, the energy 
threshold plays a key role for directional detection. It 
is worth emphasizing that it is the lowest energy at 
which both the initial direction and the energy of the 
recoiling nucleus can be retrieved, which makes it even 
more challenging for directional detection. In particular, 
this "directional threshold" is higher than the standard 
energy threshold defined as the minimal kinetic energy 
required for an event to be detected. Indeed, a low en- 
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FIG. 4: Lower bound of the 3<j discovery region at 90% CL. 
in the (m x , \og 10 (a p )) plane. Black line presents the back- 
ground free case, while the blue dashed line presents the same 
region with 10 background events per year per kg. For con- 
venience, the curves of iso-number of WIMP events are pre- 
sented (dashed lines) in the case of 5, 30 and 150 WIMP 
events. 



ergy recoil in a low pressure TPC would present a short 
track length and a large angular dispersion, implying a 
loss of the direction information. The directional energy 
threshold is thus closely related to the gas pressure, 
the target choice, the readout performance as well as 
the data analysis strategy. There are two main and 
competing effects when increasing the energy threshold 
: a reduction of the number of expected WIMP events 
and a selection of the most anisotropic WIMP-induced 
recoils. 

On figure [5l we present the evolution of E(Z) as a 
function of E t h for four different WIMP masses : m x — 
5 (black solid line), 10 (red dotted line), 20 (blue long 
dashed line) and 50 GeV/c 2 (green short dashed line). 
For the sake of comparison between the different WIMP 
masses, we have chosen cross section values in order 
to have 100 WIMP events in each case at E t h = 
keV and a constant background contamination of 10 
evts/kg/year. As discussed in [6], the fact that the 
angular recoil distribution of WIMP events is more 
anisotropic at low WIMP mass comes from the fact that 
the energy threshold selects events at higher energies 
which are the most anisotropic ones. The fact that, 
at E t h = keV, the significance is higher for lighter 
WIMPs is only due to the energy spectrum which is 
more peaked at low energy and then more different 
to the background one, assumed to be flat. Then, 
increasing the energy threshold will decrease the number 
of WIMP events more rapidly in the case of light 
WIMPs than for heavier WIMPs and also enhance the 
anisotropy related to the rejection power. However, as 
it can be seen from figure the significance is always 
decreasing with the energy threshold leading to the 
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FIG. 5: Mean significance E(Z) as a function of the energy 
threshold, for four different WIMP masses : m x = 5 (black 
solid line) , 10 (red dotted line) , 20 (blue long dashed line) and 
50 GeV/c 2 (green short dashed line). We have considered a 
constant background rate of Rb = 10 evts/kg/year and differ- 
ent values of <j v to get fi s = 100 events for each WIMP mass 
at E th = keV. 



FIG. 6: Lower bound of the 3cr discovery region at 90% CL. 
in the (m x , log 10 (cr p )) plane for two different cases: Eth = 
keV (black solid line) and Eth = 50 keV (blue long dashed 
line) without background events. For convenience, the curves 
of iso-number of WIMP events are presented in the Eth = 
keV case (black dashed lines) and in the E t h = 50 keV case 
(blue dashed lines). 



conclusion that the loss of WIMP events has a stronger 
effect on the significance than the enhancement of the 
anisotropy. Then, contrarily to what was predicted 
by [34j, when looking at the whole directional event 
distribution dr R/dE r dfl r , there is no enhancement of 
the significance at low energy threshold. Of course, 
we caution that our result presented here suffers from 
energy background modeling dependence, and that for 
other assumptions, one could find a result similar to the 
one found in [34j which as been obtained by considering 
only the angular distribution of the events. 

Figure [5] presents the lower bound of the 3ct discovery 
region at 90% C.L. in the (m x , log 10 (cr p )) for two 
different cases : E t h = keV (black solid line) and 
Eth — 50 keV (blue long dashed line). For convenience, 
the curves of iso-number of WIMP events, relevant to 
each case, are presented (dashed lines). One can first 
notice that a 30 kg. year CF4 directional detection would 
allow to achieve, in the unrealistic E t h = keV case, 
a 3<7 discovery of Dark Matter down to ~ 10~ 5 pb at 
low WIMP mass (~ 10 GeV/c 2 ) and down to - 10~ 4 
pb for high WIMP mass (~ 1 TeV/c 2 ). At high WIMP 
mass, about 4 times more WIMP events are required to 
achieve a 3cr discovery. Considering a detector with a 50 
keV energy threshold results in a loss of a factor of 5 in 
sensitivity at high WIMP mass and to no sensitivity at 
all below - 20 GeV/c 2 . 

However, as one can see from figure [SI the fact that, at 
high WIMP mass, only 15 WIMP events are required 
to get a 3er discovery at 90% C.L. for E t h — 50 keV 
while 20 are required for E t h = keV comes from the 
enhancement of the anisotropy due to the high energy 
threshold. 



As a conclusion of this study aiming at quantifying 
the impact of the energy threshold on the expected sen- 
sitivity of a given directional detection, it can simply be 
stated that the lower remains as expected the better. We 
also found that a low energy threshold is compulsory to 
be sensitive to low WIMP mass. As an illustration of this 
statement, we found that a threshold of about 50 keV 
will prevent a given directional experiment, using Fluo- 
rine target material, to be sensitive to very low WIMP 
mass ^5—10 GeV/c 2 . We conclude that the energy 
threshold remains a major experimental issue for direc- 
tional detection, noticing that both the 3D track and the 
energy must be measured down to this energy. 



C. Effect of the sense recognition efficiency 

Not only should the track be 3D-reconstructed, but 
its sense should also be retrieved from the data analysis. 
As outlined in [Hj], an asymmetry between upgoing 
and downgoing tracks is expected, due to two different 
effects. First, the angular dispersion of recoiling tracks 
should result in a spatial asymmetry as the beginning 
of the track should be more rectilinear than its end. 
Second, a charge collection asymmetry is expected as 
the dE/dx in ionization is decreasing with energy at low 
recoil energy. Hence, more primary electrons should be 
generated at the beginning of the track. 
Even though several experimental progresses have been 
done [3l|-[33, sense recognition remains a key and 
challenging experimental issue for directional detection 
of Dark Matter. In particular, it should still be demon- 
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FIG. 7: Mean significance as a function of enr for two 
different cases (N s — 100, Nb = 0) (black solid line) and 
(N s = 100, N b = 100) (blue long dashed line). A WIMP 
mass of 50 GeV/c 2 has been considered for this study. 
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FIG. 8: Lower bound of the 3cr discovery region at 90% C.L. 
in the (m x , log 10 (a p )) plane for two different cases: ejjt = 
100% (black solid line) and €ht = 0% (blue long dashed line) 
without background events. For convenience, the curves of 
iso-number of WIMP events are presented (dashed lines) in 
the case of N 3 = 5, 30 and 120. 



strated that sense recognition may be achieved at low 
recoil energy, where most WIMP events reside. 

In the following, we investigate the effect of no or par- 
tial sense recognition on the expected significance of a 
given directional detector. To do so, we define a sense 
recognition efficiency as : em — £g ~ £ff, where ea and 
ew corresponds to the fraction of good and wrong recon- 
structed events respectively. Then, if em = 0%, we are 
in the "flipping coin" scenario, i.e. the detector does not 
have any sense recognition capability. Then, the mod- 
ified directional event rate that takes into account this 
partial sense recognition S(R n ) is defined as : 

S(Rn) = ^12L S (+R n ) + l^Lsi-K), (10) 

where R n refers to the direction of the recoiling nucleus. 

Without sense recognition, the expected WIMP- 
induced distribution becomes less anisotropic and thus 
gets closer to the expected background event distribu- 
tion. This induces an obvious loss of discrimination 
power. Figure [7] presents the evolution of the mean 
significance E(Z) as a function of em for different cases 
(N s = 100, N b = 0) and (N s = 100, N b = 100) consid- 
ering a flat energy spectrum for the background. One 
can notice that, for both cases, the mean significance 
is monotonically increasing with em- It is also shown 
that this significance improvement is even stronger for 
values of €ht close to 100%. It means that having 
eHT equal to 0% or to 20% does not really modify the 
expected significance while having em equal to 80% or 
to 100% changes the mean significance by almost 20%. 
Also, it is shown on figure [7] that even in the presence 
of a sizeable background contamination (A = 0.5) the 
expected mean significance remains strong, i.e. above 



5<T. Indeed, the presence of background reduces the 
mean significance by about 30% all over the range of em- 

Figure |S] presents the lower bound of the 3<r discov- 
ery region at 90% C.L. in the (m x , log 10 (<7 p )) plane for 
two different cases: enr =100% (black solid line) and 
em =0% (blue long dashed line). For convenience, 
the curves of iso-number of WIMP events are presented 
(dashed lines) in the case of 5, 30 and 120 WIMP events. 
One can first notice that a 30 kg. year CF4 directional 
detection would allow to achieve, with full sense recog- 
nition capability, a 3cr discovery of Dark Matter down 
to - 10~ 5 pb at low WIMP mass (~ 10 GeV/c 2 ) 
and down to - 10~ 4 - 10~ 3 pb for high WIMP mass 
(-0.1-1 TeV/c 2 ). At high WIMP mass, about 6 times 
more WIMP events are required to achieve a 3(7 discov- 
ery. Considering a detector with no sense-recognition 
capability results in a loss of a factor of 4 at high WIMP 
mass and almost no effect at low WIMP mass. Note that, 
this statement depends on the background energy distri- 
bution (see sec. IIVF[) , 

This study leads us to the conclusion that sense recog- 
nition is not a major experimental issue. Indeed, a 
directional detector without sense recognition capabil- 
ity would still be able to achieve a 3<r discovery in the 
10~ 5 — 10~ 3 pb region. We emphasize that the same 
conclusion holds true when trying to set exclusion limits, 
as shown in .8]. This result is of major interest, as get- 
ting sense recognition from experimental data remains 
a very difficult task, probably the most difficult exper- 
imental challenge to be faced by current projects, and 
especially when trying to get a high sense recognition 
efficiency down to the energy threshold. 
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FIG. 9: Mean significance as a function of the angular reso- 
lution a~f for different cases (N 3 = 100, Nb = 0) (black solid 
line) and (N a = 100, N b = 0) (blue long dashed line). A 
WIMP mass of 50 GeV/c 2 has been considered for this study. 



FIG. 10: Lower bound of the 3<r discovery region at 90% C.L. 
in the {m x , log 10 (<r p )) plane for two different cases: cr 7 = 0° 
(black solid line) and <r 7 = 60° (blue long dashed line) without 
background events. For convenience, the curves of iso-number 
of WIMP events are presented (dashed lines) in the case of 5, 
30 and 100 WIMP events. 



D. Effect of angular resolution 

As far as directional detection is concerned, the esti- 
mation of the initial recoil direction is compulsory. This 
gives an intrinsic limitation of this detection strategy as 
recoil tracks in low pressure gaseous detectors would en- 
counter a rather large angular dispersion ("straggling" 
effect) by colliding with other nuclei of the gas. More- 
over, the gas properties imply a transverse and longitudi- 
nal diffusion of the primary electrons that will contribute 
to the angular resolution. 

Hence, data of upcoming directional detectors should suf- 
fer from rather large angular resolution. Having a finite 
angular resolution means that a recoil initially coming 
from the direction f(Q, r ) is reconstructed as a recoil com- 
ing from the direction f'(D,' r ) with a gaussian dispersion 
of width <t 7 according to the following distribution: 

K(fl r , il' r ) = e^ 2 / 2ff V ((2^) 2 / 3 a 7 erf(V2a 7 )) (11) 

where the angle 7 between the f(fi r ) and r'(Sl' r ) is defined 
by 

cos 7 = cos b r cos b' r cos(? r — l' r ) + sinb r smb' r , (12) 

where the angles I and b refer to the Galactic coordi- 
nates. 

The directional recoil rate is then given by the convo- 
lution product of the initial recoil rate and the angular 
ditribution as 

(13) 

The evolution of the expected mean significance 
E(Z) as a function of the angular resolution er 7 is 



shown on figure [H Two different cases are con- 
sidered (iV, = 100, Nf, = 0) (black solid line) and 
(N s = 100, N b = 100) (blue long dashed line) with a 
flat background. As one can see, for each configuration, 
degrading the angular resolution (i.e. increasing cr 7 ) 
results in a diminution of the expected mean significance. 
This comes from the fact that for large angular values 
of the angular resolution, the angular distribution of 
the expected WIMP events is less anisotropic and is 
getting closer to the angular distribution of background 
events, which is isotropic. From figure [9j we can see 
that an angular resolution of er 7 = 60° decreases the 
mean significance of about 30% in comparison to the 
case where cx 7 = 0°. 

The effect of angular resolution on the expected dis- 
covery potential of a directional detection experiment is 
shown on figure 1101 Indeed, we have represented in the 
(log 10 ((7p), m x ) plane, the 90% C.L. limit to reach a Dark 
Matter detection with a significance greater or equal to 
3ct considering a perfect angular resolution (black solid 
line) and ct 7 = 60° (blue long dashed line). For conve- 
nience, the black dashed line corresponds to iso- value of 
expected WIMP events for N s = 5, 30 and 100. As one 
can see, the effect of an angular resolution of 60° is very 
small at low WIMP mass due to the fact that at low 
WIMP mass the expected angular distribution is much 
more anisotropic than for heavy WIMPs. Then, at high 
WIMP mass, the effect of angular resolution is non negli- 
gible as we found a reduction in the sensitivity of about a 
factor of 3.3. This result suggests that directional detec- 
tion should still be very competitive even with a rather 
poor angular resolution. It may be noticed from figure [5] 
and figure [TO] that the effect of an angular resolution of 
60° with €ht = 100% is almost equivalent to not having 
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FIG. 11: Mean significance as a function of oe t for different 
cases (N s = 100, N b = 0) (solid lines) and (N s = 100, N b = 
100) (long dashed lines) according to two different WIMP 
masses: m x — 20 GeV/c 2 (black lines) and m x = 100 GeV/c 2 
(blue lines). 



the sense recognition capability with a perfect angular 
resolution. 

Then, as for the sense recognition study, we can conclude 
that angular resolution is not a major experimental issue. 
Indeed, a low angular resolution detector would still be 
able to achieve a 3<r discovery in the 10 ~ 5 — 10~ 3 pb 
region. 

E. Effect of the energy resolution 



is that it will enhance the expected number of WIMP 
events. Indeed, as most of the WIMP events lie at low 
energy, if the energy resolution (<Je t ) is large, low energy 
events could pass the energy threshold and be detected. 
The second effect is that it will smooth the energy dis- 
tribution leading to a WIMP event energy distribution 
closer to the expected background one (flat) and then re- 
duce the discrimination power. 

The effect of energy resolution is shown on figure QT] 
which presents the evolution of the mean expected sig- 
nificance E(Z) as a function of <Je t - The latter presents 
the results obtained in the cases of a 20 GeV/c 2 (black 
lines) and 100 GeV/c 2 (blue lines) WIMP mass consider- 
ing (N s = 100, N b = 0) (solid lines) and (N s = 100, N b = 
100) (long dashed lines). We have considered value for 
<je t ranging from keV, i.e. perfect energy resolution, 
to 10 keV. As the recoil energy range of the detector is 
taken, for this example, between 5 and 50 keV, oe t = 10 
keV corresponds to an extreme and unrealistic energy res- 
olution of 200% at 5 keV and 20% at 50 keV. However, as 
one can see from figure [Til even with an extremely large 
energy resolution, the mean expected significance is only 
slightly degraded compared to the case of a perfect en- 
ergy resolution, i.e. 15% in the case m x — 20 GeV/c 2 
and {N s = 100, N b = 0) (black solid line). From figure[UJ 
we can also appreciate the fact that the effect of energy 
resolution on the significance of a Dark Matter detection 
is negligible even with a 50% background contribution 
and for a 100 GeV/c 2 WIMP mass. 
As a conclusion of this study, we have shown that en- 
ergy resolution is a meaningless experimental issue for 
the Dark Matter sensitivity of a given directional detec- 
tor. 



Low pressure gaseous TPC detectors allow the ioniza- 
tion energy to be measured. It should then be converted 
to a recoil energy thanks to the knowledge of the ioniza- 
tion quenching factor. The measurement of this quantity 
is a major issue for directional detection with gaseous 

tpc ii m. 

So far, we have considered a directional detector with a 
perfect energy resolution. In this section, we present the 
effect of a non-perfect energy resolution on the expected 
sensitivity of a directional detector. As outlined above, 
most directional detection experiments are TPCs, they 
will hence be affected by energy resolutions much higher 
than the ones from cryogenic detectors for instance. 
The energy resolution on the recoil energy aE r is taken 
into account in the expected directional recoil rate dis- 
tribution using the following procedure: 



d 2 R 
dQ r dE, 
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2no(E' r ) 
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(14) 



>dE' 



Having a non-perfect energy resolution will lead to two 
different effects that are in competition. The first one 



F. Effect of the background energy model 

The last experimental issue to be discussed is the ef- 
fect of energy background modeling on the estimation of 
the sensitivity of a directional detector. In the previous 
sections, we have considered a directional event rate for 
background events flat in energy. Then, in this section, 
we consider an exponential background energy spectrum 
in the form of, 



dR 
~dK 



back 



1 x CXp (-EV/ffback) £ lg s 
-Eback exp (-Eth/Ehack) 



where -Eback refers to the slope of the background en- 
ergy distribution and if -Eback — > +oo we recover a flat 
energy spectrum. Obviously, the worst scenario is when 
the expected WIMP event energy distribution and the 
background one are the same. Indeed, in such case, 
there is no more discrimination between WIMP events 
and background events based on the energy information. 
For concreteness, we examplify this statement by consid- 
ering the case -Et>ack = 17.4 keV, corresponding to the 
slope expected for WIMP-induced recoil energy distri- 
bution when considering a WIMP mass of 50 GeV/c 2 
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FIG. 12: Mean significance as a function of the background 
slope -Eback for 5 different values of sense recognition efficien- 
cies (from top to bottom): e H r = 100%, 75%, 50%, 25% and 
0%. A WIMP mass of 50 GeV/c 2 with 100 expected WIMP 
events and no background contamination has been considered 
in this study. 



FIG. 13: Mean significance as a function of the background 
slope -Eback for 5 different values of angular resolution (from 
top to bottom): cr 7 = 0°, 15°, 30°, 45° and 60°. A WIMP 
mass of 50 GeV/c 2 with 100 expected WIMP events and no 
background contamination has been considered in this study. 



and Fluorine target. On figure [12] and [T3l one can see 
the evolution of E(Z) as a function of £back for differ- 
ent values of cht ^ n d er 7 . The common feature is that 
E(Z) is maximal at very low value of -Eb ac k corresponding 
to a very steep background energy spectrum, then de- 
creases to a minimal value located around -Eback = 17.4 
keV and then increases asymptotically to the case where 
£"back —*>■ +oo corresponding to a flat banckground en- 
ergy spectrum. An interesting point that should be high- 
lighted is the fact that the effect of -EWk on the expected 
sensitivity depends strongly on the angular performance 
of the detector. Indeed, the minimal value of E(Z) at 
-Eback = 17.4 keV is strongly degraded when decreasing 
the sense recognition efficiency or increasing the angu- 
lar resolution. This is simply explained by the fact that, 
for this particular value of the background slope exactly 
equal to the WIMP-induced event slope, the discrimina- 
tion between WIMP and background events relies only 
on the angular information. 

Finally, the effect of energy background modeling on 
the expected discovery potential of a given directional 
detector is shown on figure [T3J We have presented, on 
the (m x , log 10 (<7 p )) plane, the lower bound of the 3a 
discovery region at 90% C.L. in three different cases: 
flat background energy spectrum (-finback — ► +oo) with 
(ht = 100% (black solid line), E hack = 17 keV with 
(ht = 100% (blue long dashed line) and £"back = 17 
keV with €ht — 0% (red dotted line). Taken at face 
value, these results suggest that the effect of the back- 
ground energy distribution is negligible in the case of a 
high performance tracking device but can be significant if 
the angular resolution and/or sense recognition efficiency 
are poor. Indeed, as one can see from the red dotted line 
in figure 1141 the sensitivity is degraded by almost one 
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FIG. 14: Lower bound of the 3a discovery region at 90% 
C.L. in the (m x , log 10 (a p )) plane for three different cases: 
(-Eback +00 keV, e H T = 100%) (black solid line), (i? b ack = 
17.4 keV, £ht = 100%) (blue long dashed line) and (i?back = 
17.4 keV, eht = 0%) (red dotted line) without background 
events. For convenience, the curves of iso-number of WIMP 
events are presented (dashed lines) in the case of 5, 30 and 
220 WIMP events. 



order of magnitude when the WIMP mass corresponds 
to the considered value of -Eback and if the detector does 
not have sense recognition. 

V. PROSPECT 

In this last section, we are interested in evaluating 
the discovery potential of two very different directional 
experiments referred to as detector A and B, see Tab. HT1 
Both of them are considered as 10 kg of CF4 Time 



13 



Projection Chambers operated during three years with 
three dimensional track reconstruction. High perfor- 
mance characteritics are considered for detector A, while 
detector B stands for the pessimistic case. Indeed, the 
detector A will be considered with an energy threshold 
of 5 keV, with an angular resolution of 20°, a 100% sense 
recognition efficiency and no background contamination; 
while detector B is characterized by an energy threshold 
of 20 keV, an angular resolution of 50°, no sense recog- 
nition capability and a high background contamination 
of 10 cvcnts/kg/year uniformly distributed in energy. 
This way, detector B can be considered as a pessimistic 
scenario and detector A as an optimistic one that could 
also be interpreted as the ultimate directional detector. 





E th [keV] 


Rb [evts/kg/year] 


a j [°] 


£HT [%] 


Detector A 


5 





20 


100 


Detector B 


20 


10 


50 






TABLE II: Characteristics of the two detector A and B con- 
sidered. 



On figure [T5] we present the lower bound of the 3cr 
discovery region at 90% C.L. in the (m x ,log 10 (a p )) 
plane corresponding to the detectors A (red solid line) 
and B (blue solid line). For comparison with other direct 
searches, we have reported current limits from SIMPLE 
[3 (black solid line), COUPP [H[ (black dashed line) 
and KIMS [l6| (black dotted line). All limits are given 
in the pure-proton approximation (86|. To evaluate the 
discovery potential of Dark Matter corresponding to 
the detectors A and B we have reported the theoretical 
region obtained within the framework of the constrained 
minimal supersymmetric model, taken from [87|, as the 
green contour. 

The first result that can be inferred from figure [15] 
is that the sensitivity of detector A is about one order 
of magnitude below the one from detector B for WIMP 
masses above 50 GeV/c 2 and about three order of mag- 
nitude for m x = 10 GeV/c 2 . This highlights the need 
for detector optimization and justify a substantial exper- 
imental effort on directional detection to reach a high per- 
formance detector. It also means that, contrarily to de- 
tector A, detector B will not be sensitive to light WIMP 
below ~ 20 GeV/c 2 . The second result is that both de- 
tector A and B are competitive as they are 1 and 2 orders 
of magnitude below current exclusion limits. Moreover, 
both of them should be able to reach a large part of the 
predicted supersymmetric model which is motivated both 
by particle physics and cosmological constraints. Hence, 
even a low performance directional detector, with a 30 
kg. year exposure, could allow to identify WIMP events 
as such by recovering, with a high significance, the main 
incoming direction of the events. 




^IMPmasstGeV/c 2 ) 



FIG. 15: Lower bound of the 3<r discovery region at 90% C.L. 
in the (m x , log 10 (<r p )) plane for the two detectors A (red solid 
line) and B (blue solid line), considering a 30 kg. year expo- 
sure. The theoretical region, obtained within the framework 
of the constrained minimal supersymmetric model, is taken 
from [83] and shown as the green contour. Constraints from 
collider data and relic abundance are accounted for. Exclu- 
sion limits from direct searches are also shown: SIMPLE 13] 
(black solid line), COUPP QJ] (dotted line) and KIMS [11] 
(dashed line). 



VI. CONCLUSION 

In this paper, we have shown that the use of a pro- 
file likelihood ratio test statistic is very well suited when 
trying to estimate the sensitivity of a given directional 
experiment. Indeed, it allows us to propagate, using 
a frequentist approach, some astrophysical uncertainties 
that are very important in the field of direct searches of 
Dark Matter. This way, we have studied one by one the 
impact on the sensitivity of most of the experimental is- 
sues related to directional detection in the goal of direc- 
tional detection optimization. We have then been able 
to elaborate a weighted wish list where we found that 
the energy threshold and the background contamination 
are the most dominant effect. About the effect of angu- 
lar resolution and sense recognition efficiency, we found 
that they could affect strongly (about a factor of 4-5) the 
sensitivity of a directional experiment especially at high 
WIMP mass m x > 100 GeV/c 2 . On the other hand, we 
found that the energy resolution of the detector affects 
the sensitivity only in a negligible way, even in extreme 
and unrealistic cases. However, detector commissioning 
is compulsory to derive relevant exclusion limits or dis- 
covery regions. 

Finally, we found that directional detection with a 10 kg 
CF4 Time Projection Chamber should be able to reach an 
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important part of supersymmetric models and be com- 
petitive with current experimental limits, even in the case 
of a directional experiment with very low performance. 
Hence, we believe that directional detection of Dark Mat- 
ter is a very promising direct search for Dark Matter that 
should be able to clearly authenticate a genuine positive 
detection with a high significance. 
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